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(57) Abstract 

Methods and apparatus for determining leak and respiratory airflow are disclosed. A pressure sensor (34) and a differential pressure 
sensor (32) have connection with a pneumotach (24) to derive instantaneous mask pressure and airflow respectively. A microcontroller (38) 
estimates a non-linear conductance of any leak path occuring at a mask (12) as being the low pass tittered instantaneous airflow divided by 
the low pass filtered square root of the instantaneous pressure. The instantaneous leak flow is then the conductance multiplied by the square 
root of the instantaneous pressure, and the respiratory airflow is calculated as being the instantaneous airflow minus the instantaneous leak 
flow. The time constants for the low pass Altering performed by the microcontroller (38) can be dynamically adjusted dependent upon 
sudden changes in the instananeous leak flow. 
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Determination of Leak and Respiratory Airflow 
Field of the Invention 

The invention relates to methods and apparatus for the determination of leakage 
airflow and true respiratory airflow, particularly during mechanical ventilation. 

The airflow determination can be for a subject who is either spontaneously or 
non-spontaneously breathing, or moves between these breathing states. The invention is 
especially suitable for, but not limited to, normally conscious and spontaneously 
breathing human subjects requiring long term ventilatory assistance, particularly during 
sleep. 

Background of the Invention 

In this specification any reference to a "mask" is to be understood as including 
all forms of devices for passing breathable gas to a person's airway, including nose 
masks, nose and mouth masks, nasal prongs /pi I lows and endotracheal or tracheostomy 
tubes. 

During mechanical ventilation, breathable gas is supplied for example via a 
mask, at a pressure which is higher during inspiration and lower during expiration. It 
is useful to measure the subject's respiratory airflow during mechanical ventilation to 
assess adequacy of treatment, or to control the operation of the ventilator. 

Respiratory airflow is commonly measured with a pneumotachograph placed in 
the gas delivery path between the mask and the ventilator. Leaks between the mask and 
the subject are unavoidable. The pneumotachograph measures the sum of the 
respiratory airflow plus the flow through the- leak. If the instantaneous flow through the 
leak is known, the respiratory airflow can be calculated by subtracting the flow through 
the leak from the flow at the pneumotach. 

Known methods to correct for the flow through the leak assume (i) that the 
leak is substantially constant, and (ii) that over a sufficiently long time, inspiratory and 
expiratory respiratory airflow will cancel. If these assumptions are met, the average 
flow through the pneumotach over a sufficiently long period will equal the magnitude 
of the leak, and the true respiratory airflow may then be calculated as described. 
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The known method is only correct if the pressure at the mask is constant. If 
the mask pressure varies with time (for example, in the case of a ventilator), 
assumption (i) above will be invalid, and the calculated respiratory airflow will 
therefore be incorrect. This is shown markedly in Figs, la-lf. 

5 . 

Fig. la shows a trace of measured mask pressure in bi-level CPA P treatment 
between about 4 cm H 2 0 on expiration and 12 cm H 2 0 on inspiration. Fig. lb shows 
a trace of true respiratory airflow in synchronism with the mask pressures At time = 
21 seconds a mask leak occurs, resulting in a leakage flow from the leak that is a 
o function of the treatment pressure, as shown in Fig. 1c. The measured mask flow 
shown in Fig. Id now includes an offset due to the leak flow. The prior art method 
then determines the calculated leak flow over a number of breaths, as shown in Fig. le. 
The resulting calculated respiratory flow, as the measured flow minus the calculating 
leak flow is shown in Fig. If, having returned to the correct mean value, however is 
5 incorrectly scaled in magnitude, giving a false indication of peak positive and negative 
airflow. 

Another prior art arrangement is disclosed in European Publication 
No. 0 714 670 A2, including a calculation of a pressure-dependent leak component. 
The methodology relies on knowing precisely the occurrence of the start of an 
inspiratory event and the start of the next inspiratory event. In other words, the leak 
calculation is formed as an average over a known breath and applied to a subsequent 
breath. 

This method cannot be used if the moment of start and end of the previous 
breath are unknown. In general, it can be difficult to accurately calculate the time of 
start of a breath. This is particularly the case immediately following a sudden change 
in the leak. 

Furthermore, the method will not work in the case of a subject who is making 
no respiratory efforts, and is momentarily not being ventilated at all, for example 
during an apnea, because for the duration of the apnea there is no start or end of breath 
over which to make a calculation. 

The present invention seeks to provide a determination of leak flow and true 
respiratory airflow, accounting for the variations in flow through a leak as a function of 
pressure. 
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Summary of the Invention 

The invention discloses a method for determining instantaneous leak flow at a 
mask having a leak path during mechanical ventilation, the method comprising the steps 
5 of: 

(a) determining instantaneous airflow at the mask; 

(b) determining instantaneous pressure at the mask; 

(c) estimating non-linear conductance of said leak path as the low-pass 
filtered instantaneous airflow divided by the low-pass filtered square root of the 

io x instantaneous pressure; and 

(d) determining said instantaneous leak flow to be said conductance 
multiplied by the square root of the said instantaneous pressure. 



The invention further discloses a method for determining instantaneous 
is respiratory airflow for a subject receiving breathable gas by a mask and in the presence 
of any mask leak, the method comprising the steps of: 

(a) determining instantaneous airflow at the mask; 

(b) determining instantaneous pressure at the mask; 

(c) estimating non-linear conductance of said leak path as the low pass 
20 filtered instantaneous airflow divided by the low pass filtered square root of the 

instantaneous pressure; 

(d) determining instantaneous leak flow to be said conductance multiplied 
by the square root of the said instantaneous pressure; and 

(e) calculating the respiratory airflow as the instantaneous airflow minus 
25 the instantaneous leak flow. 



The invention yet further discloses apparatus for determining respiratory 
airflow for a subject receiving breathable gas by a mask and in the presence of any 
mask leak, the apparatus comprising: 

so transducer means located at or proximate the mask to determine instantaneous 

mask airflow and pressure; and 

processing means for estimating non-linear conductance of said leak path as the 
low pass filtered instantaneous airflow divided by the low pass filtered square root of 
die instantaneous pressure, determining instantaneous leak flow to be said conductance 

35 multiplied by the square root of the said instantaneous pressure, and calculating the 
respiratory airflow as the instantaneous airflow minus the instantaneous leak flow. 

The invention yet further discloses apparatus for providing continuous positive 
airway pressure treatment or mechanical ventilation, the apparatus comprising: 
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a turbine for the generation of a supply of breathable gas; 
a gas delivery tube having connection with the turbine; 

a mask having connection to the delivery tube to supply said breathable gas to 
a subject's airway; 

transducer means located at or proximate the mask to determine instantaneous 
mask airflow and pressure; 

processor means for estimating non-linear conductance of said leak path as the 
low pass filtered instantaneous airflow divided by the low pass filtered square root of 
the instantaneous pressure, determining instantaneous leak flow to be said conductance 
multiplied by the square root of the said instantaneous pressure, and calculating the 
respiratory airflow as the instantaneous airflow minus the instantaneous leak flow; and 

control means to control the flow generator to, in turn, control the mask 
pressure and/or mask airflow on the basis of the calculated respiratory airflow. 

The invention yet further discloses a computer program for executing the steps 
referred to above. 

In one preferred form, time constants of the low pass filtering are dynamically 
adjusted dependent upon sudden changes in the instantaneous leak flow. 

Embodiments of the invention provide advantages over the prior art. There is 
no need to know when transitions between respiratory phases occurs. The 
independence from knowledge of the subject's respiratory state has the important result 
that the leak flow calculation is accurate in apneic (i.e. no flow) instances on the part of 
25 the subject or the mechanical ventilator. 

Brief Description of the Drawings 

Embodiments of the invention will now be described with reference to the 
30 accompanying drawings in which: 

Figs, la- If shows traces of pressure and airflow from which respiratory 
airflow is calculated in accordance with a prior art method; 

Figs. 2a and b show schematic diagrams of two embodiments of ventilatory 
assistance apparatus; 

35 Fig. 3 is a block flow diagram of a method for determining instantaneous 

respiratory airflow; and 

Fig. 4 shows traces of pressure, airflow and other variables from which 
respiratory airflow is calculated; 
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Fig. 5 shows a schematic diagram of ventilatory assistance apparatus of another 
embodiment; 

Fig. 6 shows a fuzzy membership function for the calculation of the extent Aj 
to which the time t z , since the most recent positive going zero crossing of the calculated 
5 respiratory airflow is longer than the expected time T,; 

Fig. 7 shows a fuzzy membership function for the calculation of the extent B, 
to which the calculated inspiratory respiratory airflow f RESP is large positive; 

Fig. 8 shows a fuzzy membership function for the calculation of the extent A E 
to which the time t ZE since the most recent negative going zero crossing in the 
10 calculated respiratory airflow is longer than the expected time T E ; 

Fig. 9 shows a fuzzy membership function for the calculation of the extent B E 
to which the respiratory airflow f RESp is large negative; and 

Fig. 10 shows the relation between an index J and time constant t. 

is Detailed Description of Preferred Embodiments 

Fig. 2a shows mechanical ventilation apparatus 10 embodying the invention. 

The subject/patient wears a nose mask 12 of any known type. The subject 
20 equally could wear a face mask or nasal prongs/pillows, or alternatively have an 
endotracheal tube or tracheostomy tube in place. A turbine/blower 14, operated by a 
mechanically coupled electrical motor 16, receives air or breathable gas at an inlet 18 
thereof, and supplies the breathable gas at a delivery pressure to a delivery tube/hose 20 
having connection at the other end thereof with the nose mask 12. Breathable gas thus 
25 is provided to the subject's airway for the purpose of providing assisted respiration, 
with the subject's expired breath passing to atmosphere by an exhaust 22 in the delivery 
tube 20, typically located proximate to the mask 12. 

A pneumotachograph 24 is placed in the delivery tube 20 between the mask 12 
so and the exhaust 22 to provide two pressure signals, P2 and ?\ y across the 
pneumotachograph, each passed by hoses 28,30 to a differential pressure sensor 32. A 
determination of the flow of gas in the mask 12 is made the differential pressure, ?2~ 
P], resulting in a flow signal f d . The mask pressure, P 2 , also is passed to a pressure 
sensor 34 by a tapped line 36 taken from the respective hose 28, to generate a delivery 
35 pressure signal, p m , output from the pressure sensor 34. 

Both the flow signal, f<j, and the pressure signal, p m , are passed to a 
microcontroller 38 where they are sampled for subsequent signal processing, typically 
at a rate of 50 Hz. 
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The microcontroller 38 is programmed to process the flow and pressure signals 
( f d> p m) t0 produce an output control signal, y 0 , provided to an electronic motor servo- 
controller 42 that, in turn, produces a motor speed control output signal, v 0 . This 
signal is provided to the motor 16 to control the rotational speed of the turbine 14 and 
provide the desired treatment pressure, P2, at the nose mask 12. 

The motor servo-controller 42 employs a negative feedback control technique 
that compares the actual delivery pressure, in the form of the signal p m> with the 
control signal y 0 . For convenience, this control stratagem may be independent of 
operation of the microcontroller 38. 

Operation of the controlling of the microcontroller 38, so far as a calculation 
of respiratory airflow is concerned, broadly is as follows. In a sampled manner, the 
conductance of any mask leak is calculated, then the instantaneous flow through the 
leak is calculated. The flow through the leak is subtracted from the total mask flow to 
calculate the true instantaneous respiratory airflow. 

Fig. 2b shows an alternative embodiment of a system for determining true 
respiratory airflow during mechanical ventilation. The mechanical ventilation system 
10' of Fig. lb differs from that of Fig. la firstly in that the microcontroller 38 plays no 
part in control of the ventilator 50, rather only receives and data processes the 
electrically transduced mask pressure and flow signals p m , to determine and generate 
the instantaneous respiratory flow fRESP- The ventilator 50 has an internal drive signal 
provided by an oscillator 44. The motor servo controller also may or may not receive 
the mask pressure signal p m as a form of feedback control. Indeed, the ventilator 50 
can be realised by any convenient form of known generic ventilation device. 

The controlling software resident within the microcontroller 38 performs the 
following steps in determining the respiratory airflow as broadly described above, as 
also shown in the flow diagram of Fig. 3. 

The word "average" is used herein in the most general sense of the result of a 
low pass filtering step, and is not confined to an arithmetic mean. 

1. Repeatedly sample the mask airflow fj to give a sampled signal f|viASK> and the 
mask pressure P m to give a sampled signal PjviASK* for example at intervals of T=20 
milliseconds. (Steps 50,52). 
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2. Calculate the average leak. LP(L), as being the result of low pass filtering the 
airflow, f^fASK- witn * time constant of 10 seconds. (Step 54). 

3. Calculate the average of the square root of the mask pressure, LP(>/7 MASR ), as being 
5 the result of low pass filtering the square root of the mask pressure, PmaSK- with a 

time constant of 10 seconds. (Step 56). 

4. Calculate the conductance, G, of any leak (Step 58), from the equation: 
10 G = LP(L) / LP(VTW) . 

5. Calculate the instantaneous leak flow, fLEAK, through the leak (Step 60), from the 
equation: 



W = gTp 



35 



MASK 



If there is no leak flow, the value of LP(L) will be equal to zero, as will G and hence 
fLEAK- Thus the methodology is valid also where leak is equal to zero - no leak. 

20 At this juncture the leak flow has been determined, such as would be desired for a leak 
flow detector. If desired, the instantaneous respiratory airflow can be subsequently 
determined by the following step. 

6. Calculate the instantaneous respiratory airflow, fRESP, by subtracting the 
25 instantaneous leak from the mask flow (Step 62): 

fRESP = fMASK- fLEAK 

Figs. 4a-4h illustrate the methodology of the embodiment described above with 
reference to Fig. 2b. At time, t = 21 sec, ^continuing leak of approximately 1 l/sec is 
30 introduced. Fig. 4e shows the mean mask flow. Fig. 4f represents the calculated 
conductance G, from which the mask leak flow can be estimated as shown in Fig. 4g. 
Finally, Fig. 4h shows how the calculated respiratory airflow recovers within 
approximately 30 seconds, and, importantly, gives the correctly scaled (true) magnitude 
of airflow. 



With regard to setting the instantaneous output signal y 0 , the microcontroller 
broadly executes the following steps: 
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7. If the calculated true respiratory airflow fRESP is above a threshold, for example 
0.05 L/sec, y 0 is set to a value corresponding to an inspiratory pressure, PjNSP- 
Otherwise y 0 is set to a value corresponding to an expiratory pressure, PgXP- In 
general, Pj NS p is higher than P EX p, but in the case of continuous positive airways 
pressure, Pgxp may be equal to PiNSP- (Step 66). 

It is to be understood that many other methods of determining y 0 from fMASK 
may be used in step 7, for example as described in the text Principles and Practice of 
Mechanical Venti lation, edited by Martin J. Tobin (McGraw Hill Inc, 1994). 

In order to control ventilation, it is necessary to measure the subject's 
ventilation. In the presence of a leak, the ventilation delivered by the assisted 
ventilation apparatus is greater than the ventilation delivered to the subject. Known 
devices which servo-control ventilation cope with this by collecting the exhaled air 
stream with a complex system of valves, and then measuring the exhaled ventilation. 
This is inappropriate for devices for use in a domestic setting during sleep, because of 
the attendant weight, complexity, and expense. The embodiment described compensates 
for the leak by continuously measuring the nonlinear conductance of the leak, and 
allowing for the instantaneous flow through the leak as a function of pressure. 

Fig. 5 shows an alternate arrangement for ventilatory assistance apparatus 10' 
embodying the invention. In this arrangement, the pneumotachograph 24' is interposed 
between the turbine 14 and the delivery hose 20. 

This arrangement removes the pressure sensing hoses and pneumotachograph 
from the region of the mask 12. The pressure at the mask, PmaSK> is calculated from 
the delivery, pressure at the turbine 14, and from the pressure drop down the air 
delivery hose 20, which for any particular delivery hose is a known function of the flow 
at the pneumotachograph 24. Further, the microcontroller 38 must also calculate the 
flow through the mask from the flow at the turbine 14 less the flow through the exhaust 
22, which for any particular exhaust is a known function of the pressure at the mask 12. 

In more detail, this involves the steps of, firstly measuring the pressure P3 at 
the turbine 14 with the pressure sensor 34 to produce an electrical signal p t . Next the 
differential pressure P4-P3 is measured across the pneumotachograph 24' by the 
differential pressure sensor 32 to produce an electrical signal f t . In a sampled manner, 
p t and f t are digitized to yield the sampled turbine pressure and flow signals 
p TURBINE and F TURBINE- 
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The pressure at the mask PmaSK and toe sampled airflow at the mask fMASK 
12 are calculated from the turbine pressure PtURBINE and *e flow at the outlet of the 
turbine F-puRBINE & follows: 

1. Calculate the pressure drop APjuBE down the air delivery tube 20, from the 
flow at the outlet of the turbine FtURBINE : 

APTUBE = sign (F TURBINE ) x K| (F TURB[NE )2 + K 2 F TURBINE 

where K\ and K 2 are empirically determined constants, and sign (x) is 1 for 
x > 0 and -1 otherwise. 



15 



20 



2. Calculate the pressure at the mask, PmaSK> as th e pressure at the turbine 
P TURBINE les s the pressure drop AP TU BE down the air delivery tube.20: 

p MASK = p TURBINE * APTUBE- 

3. Calculate the flow, f E XHAUST> through the exhaust 22, from the pressure at 
the mask P\4ASK : 



f EXHAUST = si g n ( p MASK) x ^3 v/abs(P MASK ) 
where K3 is determined empirically. 

25 4. Calculate the flow, fMASK* int0 the mas k 12 as the flow at the turbine 14 less 
the flow through the exhaust 22: 



f MASK = ^TURBINE - ^EXHAUST- 

The foregoing embodiments describe low-pass filtering of both the 
instantaneous airflow and the square root of the instantaneous pressure with a time 
constant x of 10 seconds. This time constant, t, can be advantageously dynamically 
adjustable. 



If the conductance of the leak suddenly changes, then the calculated 
conductance will initially be incorrect, and will gradually approach the correct value at 
a rate which will be slow if the time constant of the low pass filters is long, and fast if 
the time constant is short. Conversely, if the impedance of the leak is steady, the 
longer the time constant the more accurate the calculation of the instantaneous leak. 
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Therefore, it is desirable to lengthen the time constant if it is certain that the leak is 
steady, reduce the time constant if it is certain that the leak has suddenly changed, and 
to use intermediately longer or shorter time constants if it is intermediately certain that 
the leak is steady. 

If there is a large and sudden increase in the conductance of the leak, then the 
calculated respiratory airflow will be incorrect. In particular during apparent 
inspiration, the calculated respiratory airflow will be large positive for a time that is 
large compared with the expected duration of a normal inspiration. Conversely, if there 
is a sudden decrease in conductance of the leak, then during apparent expiration the 
calculated respiratory airflow will be large negative for a time that is large compared 
with the duration of normal expiration. 

Therefore, an index of the degree of certainty that the leak has suddenly 
changed is derived, such that the longer the airflow has been away from zero, and by a 
larger amount, the larger the index; and the time constant for the low pass filters is 
adjusted to vary inversely with the index. In operation, if there is a sudden and large 
change in the leak, the index will be large, and the time constant for the calculation of 
the conductance of the leak will be small, allowing rapid convergence on the new value 
of the leakage conductance. Conversely, if the leak is steady for a long time, the index 
will be small, and the time constant for calculation of the leakage conductance will be 
large, enabling accurate calculation of the instantaneous respiratory airflow. In the 
spectrum of intermediate situations, where the calculated instantaneous respiratory 
airflow is larger and for longer periods, the index will be progressively larger, and the 
time constant for the calculation of the leak will progressively reduce. For example, at 
a moment in time where it is uncertain whether the leak is in fact constant, and the 
subject merely commenced a large sigh, or whether in fact there has been a sudden 
increase in the leak, the index will be of an intermediate value, and the time constant 
for calculation of the impedance of the leak will also be of an intermediate value. One 
advantage is that some corrective action will occur very early. 

Another advantage is that there is never a moment where the leak correction 
algorithm is "out of control" and needs to be restarted, as described for prior art 
European Patent Publication No. 0 714 670 A2. 

In a preferred embodiment, the above index is derived using fuzzy logic. The 
fuzzy extent Aj to which the airflow has been positive for longer than expected is 
calculated from the time l Z \ sin ce the last positive-going zero crossing of the calculated 
respiratory airflow signal, and the expected duration Tj of a normal inspiration for the 
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particular subject, using the fuzzy membership function shown in Fig. 6. The fuzzy 
extent Bj to which the airflow is large and positive is calculated from the instantaneous 
respiratory airflow using the fuzzy membership function shown in Fig. 7. The 
instantaneous index I { of the degree of certainty that the leak has suddenly increased is 
5 calculated as the fuzzy intersection (lesser) of A| and Bj. 

Comparable calculations are performed for expiration as follows. The fuzzy 
extent A E to which the airflow has been negative for longer than expected is calculated 
from the time t ZE since the last negative-going zero crossing of the calculated 

io respiratory airflow signal, and T E , the expected duration of a typical expiration for the 
particular subject, using the membership function shown in Fig. 8. The fuzzy extent 
B E to which the airflow is large negative is calculated from the instantaneous 
respiratory airflow using the fuzzy membership function shown in Fig. 9. The 
instantaneous index 1 E of the degree of certainty that the leak has suddenly decreased is 

is calculated as the fuzzy intersection of A E and B E . 

The instantaneous index I of the extent to which there has been a sudden 
change in the leak (either an increase or a decrease) is calculated as the fuzzy union 
(larger) of indices I j and 1 E . The instantaneous index I is then passed through a peak 

20 detector followed by a low pass filter with a time constant of, for example 2 seconds, to 
yield the desired index J. Thus if index I becomes momentarily large, index J will be 
initially large and remain so for a few seconds, The time constant x for the low pass 
filters used in the calculation of the conductance of the leak is then adjusted to vary 
inversely with the index J, as shown in Fig. 10. For example, if the expected duration 

25 of a normal respiratory cycle were 4 seconds the time constant is set to 10 seconds if 
the index J is zero, (corresponding to complete certainty that the leak is steady), and to 
1 second if the index J is unity (corresponding to complete certainty that the leak is 
suddenly changing), and to intermediate values for intermediate cases. 

30 The embodiments described refer to apparatus for the provision of ventilatory 

assistance, however, it is to be understood that the invention is applicable to all forms 
of mechanical ventilation and apparatus for the provision of continuous positive airway 
pressure treatment. The apparatus can be for the provision of a constant treatment 
pressure, multi-level (IPAP and EPAP) treatment or autosetting (adjusting) treatment or 

as other forms of mechanical ventilation, including Proportional Assist Ventilation (PAV) 
as taught by M Younes in the above-noted text. 
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The methodology described can be implemented in the form of a computer 
program that is executed by the microcontroller described, or by discrete combinational 
logic elements, or by analog hardware. 
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1. A method for determining instantaneous leak flow at a mask having a 
leak path during mechanical ventilation, the method comprising the steps of: 
5 (a) determining instantaneous airflow at the mask; 

(b) determining instantaneous pressure at the mask; 

(c) estimating non-linear conductance of said leak path as the low-pass 
filtered instantaneous airflow divided by the low-pass filtered square root of the 
instantaneous pressure; and 

10 (d) determining said instantaneous leak flow to be said conductance 

multiplied by the square root of the said instantaneous pressure. 

2. A method for determining instantaneous respiratory airflow for a 
subject receiving breathable gas by a mask and in the presence of any mask leak, the 
is method comprising the steps of: 

(a) determining instantaneous airflow at the mask; 

(b) determining instantaneous pressure at the mask; 

(c) estimating non-linear conductance of said leak path as the low pass 
filtered instantaneous airflow divided by the low pass filtered square root of the 

20 instantaneous pressure; 

(d) determining instantaneous leak flow to be said conductance multiplied 
by the square root of the said instantaneous pressure; and 

(e) calculating the respiratory airflow as the instantaneous airflow minus 
the instantaneous leak flow. 



3. A method as claimed in claim 2, whereby the time constants for said 
low pass filtering are dynamically adjustable dependent upon sudden changes in said 
instantaneous leak flow. 

4. A method as claimed in claim 3, whereby said dynamic adjustment 
comprises the further steps of: 

deriving an index of the extent to which said conductance has changed 
suddenly; and 

changing said time constants in an opposite sense to a corresponding change in 
35 said index. 

5. A method as claimed in claim 4, whereby said index is derived by the 

steps of: 
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from said calculated respiratory airflow, determining the extent to which the 
absolute magnitude of calculated airflow is larger than expected for longer than 
expected. 

6. A method as claimed in claim 2, whereby steps (a) and (b) comprise: 
measuring airflow and pressure in a gas delivery circuit coupled to said mask; 
calculating the pressure drop along the delivery circuit to the mask as a 
function of said delivery circuit airflow; and 

calculating a derived said instantaneous mask pressure as the measured delivery 
circuit pressure less the pressure drop; and 

calculating the airflow through an exhaust of the mask as a function of the 
derived instantaneous mask pressure; and 

calculating a derived said mask airflow as the measured delivery circuit airflow 
minus the exhaust airflow. 



7. Apparatus for determining respiratory airflow for a subject receiving 
breathable gas by a mask and in the presence of any mask leak, the apparatus 
comprising: 

transducer means located at or proximate the mask to determine instantaneous 
mask airflow and pressure; and 

processing means for estimating non-linear conductance of said leak path as the 
low pass filtered instantaneous airflow divided by the low pass filtered square root of 
the instantaneous pressure, determining instantaneous leak flow to be said conductance 
multiplied by the square root of the said instantaneous pressure, and calculating the 
respiratory airflow as the instantaneous airflow minus the instantaneous leak flow. 

8. Apparatus as claimed in claim 7, wherein the time constants for said 
low pass filtering are dynamically adjustable dependent upon sudden changes in said 
instantaneous leak flow. 

9. Apparatus as claimed in claim 8, wherein said processor means 
dynamically adjusts the time constants by deriving an index of the extent to which said 
conductance has changed suddenly, and changing said time constants in an opposite 
sense to a corresponding change in said index. 

10. Apparatus as claimed in claim 9, wherein said processor means 
derives said index from said calculated respiratory airflow by determining the extent to 
which the absolute magnitude of calculated airflow is larger than expected for longer 
than expected. 
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11. Apparatus as claimed in claim 7, wherein said transducer means 
comprises a pneumotachograph coupled to a differential pressure transducer. 

5 12. Apparatus as claimed in claim 11, wherein said pneumotachograph is 

located between the mask and the mask exhaust. 

13. Apparatus as claimed in claim 11, wherein said transducer means is 
located in a gas delivery circuit connected with said mask and remote from said mask. 

10 

14. Apparatus for providing continuous positive airway pressure treatment 
or mechanical ventilation, the apparatus comprising: 

a turbine for the generation of a supply of breathable gas; 
a gas delivery tube having connection with the turbine; 
15 a mask having connection to the delivery tube to supply said breathable gas to 

a subject's airway; 

transducer means located at or proximate the mask to determine instantaneous 
mask airflow and pressure; 

processor means for estimating non-linear conductance of said leak path as the 
20 low pass filtered instantaneous airflow divided by the low pass filtered square root of 
the instantaneous pressure, determining instantaneous leak flow to be said conductance 
multiplied by the square root of the said instantaneous pressure, and calculating the 
respiratory airflow as the instantaneous airflow minus the instantaneous leak flow; and 
control means to control the flow generator to, in turn, control the mask 
25 pressure and/or mask airflow on the basis of the calculated respiratory airflow. 

15. Apparatus as claimed in claim 14, wherein the time constants for said 
low pass filtering are dynamically adjustable dependent upon sudden changes in said 
instantaneous leak flow. 

30 

16. Apparatus as claimed in claim 15, wherein said processor means 
dynamically adjusts the time constants by deriving an index of the extent to which said 
conductance has changed suddenly, and changes said time constants in an opposite sense 
to a corresponding change in said index. 

35 

17. Apparatus as claimed in claim 16, wherein said processor means 
derives said index from said calculated respiratory airflow by determining the extent to 
which the absolute magnitude of calculated airflow is larger than expected for longer 
than expected. 
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18. A computer program for determining instantaneous respiratory airflow 
for a subject receiving breathable gas by a mask and in the presence of any mask leak, 
the program receiving input data of instantaneous airflow and pressure at the mask, and 
comprising the computational steps of: 

(a) determining instantaneous airflow at the mask; 

(b) determining instantaneous pressure at the mask; 

(c) estimating non-linear conductance of said leak path as the low pass 
filtered instantaneous airflow divided by the low pass filtered square root of the 
instantaneous pressure; 

(d) determining instantaneous leak flow to be said conductance multiplied 
by the square root of the said instantaneous pressure; and 

(e) calculating the respiratory airflow as the instantaneous airflow minus 
die instantaneous leak flow. 
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